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Abstract 

Two recently proposed concepts to improve the perturbative calculation of 
exclusive amplitudes, gluonic radiative corrections (Sudakov factor) and con- 
finement size effects (intrinsic transverse momentum) are combined to study 
the neutron magnetic form factor in the space-like region. We find that nu- 
cleon distribution amplitudes modelled on the basis of current QCD sum rules 
indicate overlap with the existing data at the highest measured values of mo- 
mentum transfer. However, sizeable higher-order perturbative corrections (K- 
factor) and/or higher-twist contributions cannot be excluded, although they 
may be weaker than in the proton case. 



In a recent paper |]l| we have studied the space-hke proton form factor within a theo- 
retical scheme proposed by Li and Sterman [^, which takes into account gluonic radiative 
corrections in the form of a Sudakov factor . This scheme naturally generalizes the stan- 
dard hard scattering picture (HSP) of Brodsky and Lepage |Q — commonly used to calculate 
exclusive reactions within perturbative QCD — by taking into account the transverse mo- 
mentum of the partons. 

A major point in our proton form-factor analysis was to show that proper treatment 
of the a^-singularities demands the imposition of an appropriate infrared (IR) cut-off to 
render the form-factor calculation both finite and insensitive to the inclusion of the soft 
region of phase space. This is in contrast to the pion case [Q, where a "natural" IR cut-off 
appears in the form of the interquark separation. Considering in detail optional IR cut-off 
prescriptions [§-0]; we found that maximum IR protection is provided by introducing as a 
common IR cut-off in the Sudakov (suppression) factor the maximum interquark separation 
("MAX" prescription |1|]). The underlying physical idea is the following: One expects that 
because of the color neutrality of a hadron, its quark distribution cannot be resolved by 
gluons with a wavelength much larger than a characteristic interquark separation scale bi. 
Thus, gluons with wavelengths large compared to the (transverse) hadron size probe the 
hadron as a whole, i.e., in a color-singlet state and decouple. As a result, quarks in such 
configurations act coherently and therefore (soft) gluon radiation is dynamically inhibited. 

The "MAX" presription not only suffices to suppress the Os-singularities, but also pre- 
serves the finiteness of the integrand of the expression for the form factor, even when 
renormalization-group (RG) evolution of the wave function is included. An additional bonus 
of the "MAX" prescription is that the proton form factor saturates, i.e., becomes insensitive 
to the contributions from large transverse separations. Admittedly, little confidence is put 
in perturbative treatments of the large-distance region, so that saturation of the form factor 
at transverse distances as low as possible is a prerequisite for a self-consistent perturbative 
calculation. 

In [H] we have pointed out that the recent numerical analysis by Li [^ of the proton 



form factor has serious drawbacks: (i) The cancellation of the as-singularities in the region 
bi-^QCD — 1 {xi fixed) is incomplete for different 6/, amounting to uncompensated singularities 
of the form 

^\n(^J—] . (1) 

\bi Aqcd J 

(ii) There is no saturation, meaning that the main form-factor contributions are accumulated 
in the "forbidden" soft region. 

A second element of our approach in |l| was the incorporation of the intrinsic transverse 
momentum in the proton wave function, following a previous work by two of us P] on the 
pion form factor. The intrinsic transverse momentum reflects confinement-size effects 
and improves the saturation behavior of the form factor. As a consequence of these ef- 
fects (Sudakov factor and intrinsic transverse momentum), the self-consistent perturbative 
contribution to the proton form factor turns out to be reduced by at least a factor of two 



compared to the existing experimental data (see, e.g., |TO|,|Tl|). This is true for a variety 
of nucleon distribution amplitudes (DA), recently determined by two of us |T2|,|TB[ on the 
theoretical basis of QCD sum rules |T^JT^. 

In the present work we extend this type of analysis to the neutron magnetic form fac- 
tor. While in our previous work the focus was on the self-consistent implementation of the 
Sudakov factor and the proper identification and matching of the scales involved, the sub- 
ject of the current effort will be on the phenomenological side. In particular, we consider 
observables involving the ratio of the neutron to the proton form factor G^.j/G%j. It has 
been discussed in [|16| (and previous references cited therein) and more recently in |jl^ that 



proposed model distribution amplitudes for the nucleon can be classified according to this 
ratio (an observable quantity) and the theoretical parameters i?4 (projection coefficient on 
the corresponding eigenfunction of the nucleon evolution equation) and the "hybridity" an- 



gle '(9, as detailed in [O. One of the crucial questions posed in the present work is whether 



the emerging pattern of solutions to the sum rules, found within the standard HSP 0|T^ to 
constitute a smooth and finite "orbit" in the (-B4, —G\i/G\,^) plane, pertains to the inclusion 



of transverse-momentum contributions. 

The starting point of our analysis is to consider the neutron magnetic form factor within 
the modified HSP: 

G^(Q') = yI {dx][dx'] J j^^j^^Y.TA^^^'^b,Q,^^)Yp{x,x',b,f^F)e-'^, (2) 

with [dx] = dxidx2dx36{l — I] Xj); Xi being the momentum fractions carried by the valence 
quarks. The Fourier-transformed hard scattering amplitudes are given by 



Ti = - Cf as{tn)as{ti2)Ko ( \/{l - xi){l - x[) Qh ) K^ ( Va^24 Q&2 ) , (3) 



Ta = - Cf a,(t2i)as(t22)^o N xix'i QhA Kq f ^xa^ Qfoz j , (4) 

where the Kq are modified Bessel functions of order and bi denotes the length of the 
corresponding transverse-distance vector. 

The renormalization scale is chosen in such a way that each hard gluon refers to its own 
individual momentum scale tji to be used in the argument of the corresponding as- The tji is 
defined as the maximum scale of either the longitudinal momentum or the inverse transverse 
separation, associated with each of the gluons: viz. 



tii= max 
^21= max 



^{l-x,){l-x[)Q,l/b, 



\Jxix'iQ,l/bi 
ti2= t22 = max 



X2X2Q,l/b2 



(5) 



Since the hard scattering amplitudes depend only on the differences of initial and final 
state transverse momenta, there are only two transverse separation vectors, namely those 
between quarks 1 and 3 and between quarks 2 and 3: 6i (= b[), 62 (= b'2). Accordingly, the 
transverse separation between quark 1 and quark 2 is 63 = 62 — ^i- 

The soft part of the form factor is given by the following expressions which contain 
linear combinations of products of the initial and final state wave functions in the transverse 
configuration space, weighted by Xj-dependent factors arising from the fermion propagators: 



^" = (l_^^)^(l_^.) { -2^^23^123 - 2^1:^2^132 + ^231^231 + ^321^321 

- ^231^132 - ^^32^231 - ^321^123 " ^m^32l} (6) 

^2" = (l_^^)(l_3./) {^231^231 + ^231^132 + ^t32^23l} 

+ 71 Vl M (2^321^321 - ^^23^123 + ^321^123 + ^t23^32l} • (7) 

[i — X3)[i — Xi) 

The Fourier transform of the wave function reads 

^123(X, 6, /ii.) = — ^==/N(/iF)$123(a;,/iF)^123(a;,6), (8) 

wherein its intrinsic /c^-dependence is parametrized according to the Gaussian 

^123(3;, b) = (47r)^exp | -— XiX^bl + X2X^hl + X1X263 \ • (9) 

We have used the convenient short-hand notation \l/i23(a;, h) = \l/(xi, 6; X2, b; X3, b) denoting 
in $i23(x, /ii?) the factorization scale of short-and large-distance contributions by hf- 

The exponentials e""^^ in (^ are the Sudakov factors responsible for the effects of gluonic 
radiative corrections. They have been calculated by Botts and Sterman [§] using resumma- 
tion techniques in the context of the renormalization group (RG) and having recourse to 
previous extensive work by Collins, Soper, and Sterman [|l^. The explicit expressions for 
the Sudakov exponents are given in 0,||. 

The analytical and numerical evaluation of the neutron form factor is performed under 
the imposition of the "MAX" IR-prescription [0 on the Sudakov factor, i.e., setting 

b = max{6i, 62, ^3} = 61 = &2 = ^3- (10) 



The results of this calculation are typified by the curves shown in Fig. |1] for the COZ DA |T^ 
(solid line), including also the intrinsic transverse momentum in two different ways: (i) by 
normalizing the probability P^q for finding three valence quarks in the neutron to unity 
(dashed line), which results to {k\Y^'^ = 271 MeV; and (ii) by setting the value of the r.m.s. 
transverse momentum equal to 600 MeV (dotted line), which implies P-^q = 0.042. Here 
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and below we throughout use the values Aqcd = 180 MeV and |/Ar| = (5.0 ± 0.3) x 10 ^ 
GeV^ |T^, the latter being the value of the nucleon DA at the origin. The momentum 



evolution of the DA is RG-controlled — provided the model DA is satisfying the nucleon 
evolution equation |^, which is true for all DAs we consider in this work. Then the nucleon 
DA can be expanded in terms of the eigenf unctions of the one-gluon exchange kernel to read 



$i23(a;,/i) = $^^3(x)^5„ 



' Osijjy 



7n//3o 



^123 (^), 



(11) 



where the notations of |T^ are adopted and $ 



123 



X 



120x1X2X3 is the asymptotic DA. 



The exponents 7^ are related to the anomalous dimensions of trilinear quark operators with 



isospin 1/2 (see ||T9l) and resemble the 6„ in the Brodsky-Lepage notation [If. Because they 
are positive fractional numbers increasing with n, higher-order terms in (pA[) are gradually 
suppressed. The constants % are given in |jl|,|^; /?o = ll — 2nf/3 = 9 for three flavors. Under 
these conditions we may use QCD sum-rule results on the moments 



Jo 



(12) 



to constrain the first few expansion coefficients -B„ (for more details see [|1^,^,|I^). 

As can be seen from Fig. |I] and comparison with |]T^, the magnitude of the neutron 
magnetic form factor with Sudakov correction is reduced by more than a factor of 2.5 
with respect to the standard HSP. This reduction is enhanced when the intrinsic transverse 



momentum is included and is pertinent to all nucleon DAs [|T2| , p!3| modelled on the basis of 
existing QCD sum-rules [0,0 (shaded band in Fig. ^. The upper region of the band is 
characterized by COZ-like ^M DAs, whereas its lower part is associated with the recently 



proposed "heterotic" DA |2^. We note that the perturbative contribution becomes self- 
consistent for momentum transfers larger than 8 GeV^ (using {k'j_)^^'^ = 271 MeV) in the 
sense that at least 50 % of the result is accumulated in regions where a^ is smaller than 0.5. 
It is important to emphasize that the neutron magnetic form factor is the first process 
calculated within the modified HSP that yields predictions which indicate overlap with the 
existing data [^, as can be seen from Fig. |^. This tentative agreement occurs at data 



points corresponding to the largest momentum transfers measured, where, incidentally, our 
theoretical calculations become self-consistent. Therefore, measurements of the neutron 
magnetic form factor beyond 10 GeV^ are extremely important in order to check the validity 
of the theoretical predictions in a more quantitative way. 

One place to test these results is in the data for the differential cross sections for elastic 
electron-proton and electron-neutron scattering CTp and (Jn, respectively. For small scattering 
angles, where the terms oc tan^(0/2) can be neglected, and for large Q^, the ratio o"n/o"p 
becomes in a slightly mo del- dependent way proportional to the square of the ratio of the 
neutron to the proton magnetic form factor. 

Combining our calculations for the proton [|l| with those presented here for the neutron, 
we can extract theoretical predictions for a^Ja^ by inputting the same set of model DAs for 
the nucleon |[l^|l^ as before. The results are shown in Fig. |^ (shaded area) in comparison 



with available data p2|. From this figure we see that the measured values of (Tn/cTp enter the 



estimated range already at Q^ ~ 8 GeV^. [The corresponding values of the ratio —G\JG\j, 
allowed by our analysis, range between -0.2 and O.5.] The fair agreement between theoretical 
predictions and data is partly deceptive owing to the fact that the self-consistent calculation 
of the leading-order perturbative contribution to the proton magnetic form factor within the 
modified HSP yields a rather small value [|ll. This missing part of the proton form factor 
can arise from many sources, e.g, from a large K-factor and/or higher-twist contributions 
and would certainly affect the width of the strip. Such contributions are also conceivable 
for the neutron form factor. 

In any case it is remarkable that the collective pattern of solutions to the QCD sum 



rules [14,15|, found within the standard HSP |12,TJ, pertains to the inclusion of transverse- 



momentum contributions comprising the Sudakov factor and those due to the intrinsic trans- 
verse momentum (see Fig. ^. Indeed, the solutions arrange themselves across an "orbit" 
in the (i?4, — G^f /G^f) plane which is somewhat shifted compared to the original one. In 
contrast to the standard HSP version, within the present context, the "orbit" is slightly 
Q^-dependent, as shown in Fig. ^ The new "orbit" at Q^ = 30 GeV^ can be characterized 
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by the empirical relation -GlJGlj = 0.426-9.91 x 10-^54 -4.27 x 10-^51 + 4.59 x IQ-^Sf, 
which complies with that found in |12[ . The dashed line in Fig. ^ represents a similar fit for 



Q2 _ j^gs QgY2_ ^^Q observe that with increasing momentum transfer, the "orbit" within 
the modified HSP transmutes into that of the standard HSP. We note that the coefficient B4 
projects onto the eigenfunction $4(xj) and hence provides an effective measure to account for 
the antisymmetric content of the nucleon DA, since the other antisymmetric eigenfunctions 



are offset by this term [|12 



In summary, in this letter the modified HSP has been applied to the neutron magnetic 
form factor for the first time. In contrast to other cases (e.g., pion and proton electromagnetic 
form factors), the band of predictions obtained with the set of model DAs for the nucleon 
indicates overlap with the experimental data at the largest measured values of momentum 
transfer, where the theoretical predictions become self-consistent. Nevertheless, it is likely 
that in order to improve agreement with the data, several additional contributions have to be 
included: Perturbative higher-order corrections may give rise to a rather large K-factor of the 
order of 2 multiplying the leading-order contribution, as found for other large-momentum 
transfer processes |23|. However, for the case of the pion form factor, already existing 
calculations P^p5[ of the K-factor to one-loop order indicate that with an appropriate 



choice of the renormalization point, the actual value of the K-factor is rather small, i.e., 
of order unity. On the other hand, still unestimated contributions due to higher twists are 
presumably sizeable in the experimentally accessible region and may also be important. 
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FIGURES 
FIG. 1. The influence of the intrinsic transverse momentum on the neutron magnetic form 



factor within the modified HSP. The curves shown are obtained for the COZ |14| DA by impos- 
ing the "MAX" prescription including evolution. The solid line represents the results without 
A;_i_-dependence, whereas the dashed and dotted lines are obtained with (/c^)^'^ = 271 MeV and 
600 MeV, respectively. 

FIG. 2. The neutron magnetic form factor vs. Q^ . The theoretical results are obtained using 
the "MAX" prescription including evolution and normalizing the wave functions to unity. The 
shadowed strip indicates the range of predictions derived from the set of DAs determined in |l2| in 
the context of QCD sum rules (see text). The solid (dashed, dotted) line corresponds to the COZ 
(heterotic, optimized COZ) DA (cf. Fig. 4). The data are taken from pT| . 

FIG. 3. The ratio a^ /a^ of the differential elastic electron-neutron to electron-proton cross 
section vs. Q^ at scattering angles of 10°. The shaded area and model DAs for the nucleon 
correspond to those shown in Fig. 2. The data are taken from p2[ . 

FIG. 4. Relation between the ratio R of the magnetic nucleon form factors and the expansion 
coefficient B/^ of the Appell polynomial decomposition of the nucleon DA, including the effect of the 
Sudakov factor. The results are obtained at Q^ = 30 GeV^, employing the "MAX" prescription. 
The superimposed solid line is an empirical polynomial fit similar to the original one given in jl2| . 
The dashed line serves to illustrate the dependence on the momentum scale (Q^ = 10^ GeV^). 
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http://arXiv.org/ps/hep-ph/9407250v2 



This figure "figl-2.png" is available in "png" format from: 



http://arXiv.org/ps/hep-ph/9407250v2 



This figure "figl-3.png" is available in "png" format from: 



http://arXiv.org/ps/hep-ph/9407250v2 
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